Twisted speckle entities inside wavefront reversal mirrors. 
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The previously unknown property of the optical speckle pattern reported. The interference of 
a speckle with the counter-propagating phase-conjugated speckle wave produces a randomly dis- 
tributed ensemble of a twisted entities (ropes) surrounding optical vortex lines. These entities 
appear in a wide range of a randomly chosen speckle parameters inside the phase-conjugating mir- 
rors regardless to an internal physical mechanism of the wavefront reversal. These numerically 
generated interference patterns are relevant to the Brillouin PC-mirrors and to a four-wave mixing 
PC-mirrors based upon laser trapped ultracold atomic cloud. 
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PACS numbers: 42.30.Ms 42.50.Tx 42.65.Hw 42.65.Es 

Phase singularities of the optical beams attracted a 
substantial interest in a recent decades from the point 
of view of optical information processing Helical 
wavefronts had been shown to affect a processes of the 
second harmonic generation image processing with 
a photorefractive mirrors [2j, phase-conjugated reflec- 
tion via nondegenerate four-wave mixing in a cold atoms 
cloud 0]. The photons with helical wavefunction had 
been shown to possess quantized angular momentum h 
per photon 0]. The optical nonlinearities are capable 
to transfer the angular momentum from photons to an 
ensemble of ultracold atoms ||. This effect is consid- 
ered as a possible tool for the light storage in addition 
to the slow light technique based on electromagnetically 
induced transparency (E1T) [f| Recently the angu- 
lar momentum transfer to BEC cloud of sodium atoms 
via stimulated Raman scheme had been observed 0]. Of 
special interest is the nondegenerate four-wave mixing 
a cold atomic cloud and phase conjugated (PC) reflec- 
tion of the optical vortices (OV), i.e. Laguerre-Gaussian 
beams (LG) with helical wavefronts The reflection 
of the phase-conjugated LG photons from 10 7 cesium 
atoms cooled down to T « 10~ 3 K and reversal of the 
orbital angular momentum (OAM) had been interpreted 
as a consequence of internal macroscopic rotations inside 
atomic cloud |8|. Afterwards these experimental results 
were analyzed from a point of view of the angular mo- 
mentum conservation for the incident and reflected pho- 
tons from a Brillouin PC - mirror [9(. The goal of the 
current communication is to study the spiral interference 
patterns with the period of A/2 around the nodes of op- 
tical speckle pattern. These coarse interference patterns 
turn in rotation the Brillouin medium [loj producing the 
acoustical vortices carrying OAM [ill ]. 

The traditional experimental and numerical technique 
for the visualization of the optical phase singularities is 
based on mixing of a wide quasi-plane wave reference op- 
tical beam with a speckle signal [12| • The dark spots (ze- 
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ros of amplitudes) appearing in intensity distribution are 
collocated with the helical ramps in phase distribution. 
The current communication reports a theoretical descrip- 
tion of the alternative visualization technique which uses 
the interference of the speckle wave with the counter- 
propagating phase-conjugated one. In fact this technique 
is automatically implemented in a wavefront reversing 
mirrors where a phase-conjugated wave produces near 
the bright spot the Bragg grating of dielectric permit- 
tivity of the form cos((k p + k s )z). In the vicinity of the 
dark line the modulation of the light intensity is more 
complex: the Bragg grating is transformed into spirals 
of the form cos((k p + k s )z ± 2 £cf>), where (j> is the local 
azymuthal angle, £ is topological charge of the phase sin- 
gularity, k p ,k s are wavenumbers of incident (pump) and 
reflected (Stockes) waves respectively 0, In the real 
PC-mirrors the contrast of such gratings may be reduced 
due to the interference with a non-phase-conjugated com- 
ponent of the radiation. 

The novel feature compared to the previous findings 
[l2l [l3j is that the ideal (or close to an ideal) PC mirror 
visualizes in the speckle patterns a peculiar optical enti- 
ties, which might be called the ropes. Typically each such 
a rope is composed of a 2-3 optical vortex lines (figd]). In 
contrast to the well studied straight optical vortex lines 
e.g. LG [l|, [J] , the OV in a speckle field are self- twisted 
as it seen from figUl Having in mind the well known fact 
[l2T ] , that OV appear and annihilate as a pairs of whirls 
with opposite circulations 0, [T3| it is easy to conclude 
that at least two adjacent OV with opposite topological 
charges are needed to produce a rope (figE}. The mean 
length of each OV in Z direction is the Rayleigh range 
of a speckle pattern Lr w D 2 /A, where D is an average 
transverse size (in the plane X, Y) of the OV core, the 
wavelength A = 2Tr/k p . The figQ] shows the numerically 
generated fragment of a speckle pattern which contains 
three ropes each composed of a set of a vortex lines. 
Let us describe the numerical procedure for generating 
the interference pattern of the two counter-propagating 
speckle fields. 

The standard model of the phase conjugation via stim- 
ulated Brillouin scattering (SBS) is the Bragg reflection 
from the sound grating with period A/2 moving with 
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FIG. 1: Gray scale intensity plot for a small volume of 
a speckle pattern inside PC-mirror in (X, Z) -plane (Y — 
23pm). The mean transverse size of a speckle is D ~ 8pm. 
The mean longitudinal length of a speckle entity is of the order 
of the Rayleigh range Lr « D 2 /A for A ~ 1pm. The size of 
the volume is 32 pm in transverse (X, Y) directions and 128 
pm in longitudinal (Z) direction. The step of longitudinal 
modulation cos((k p + k a )z ± 2</>) induced by the interference 
of counter-propagating pump and Stockes waves is enlarged 
here 8 times for visualization purposes. The characteristic tt- 
phase-jump is clearly visible in between adjacent Bragg's cos 
-modulated roll patterns. 



the speed of sound v ac 
of momentum p = hk p 



[13j |. Due to the conservation 
= hk s + hk ac the length of the 
wavevector of sound k ac is equal to doubled length of 
wavevector of pump light k p with an accuracy about 10 -5 
[ijj]. The Doppler effect defines the small (of the order 
of 1CP 5 ) frequency shift of the reflected (Stockes) wave 
u> ac = 2u! p nv ac / c, where n is refractive index of medium, 
c is the speed of light in vacuum, uj p , lu s is the frequency 



of the pump and Stockes light respectively. 

The evolution in space of the two counter-propagating 
paraxial laser beams inside Brillouin active medium 
obeys to Maxwell equations with the cubic nonlinear po- 
larization induced by electrostrictive effect [13| . The lin- 
early polarized pump field E p moves in a positive direc- 
tion of Z — axis, the reflected Stockes field E s with the 
same polarization propagates in opposite direction(fi g.l) . 
The acoustic field Q is excited via electrostriction [la,lla|. 
The three- wave equations for the interaction of E p , E s 
with adiabatic elimination of the Q sa E p E s * are: 



dE p (z,x,y,t) 
dz 
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where 7 = p{de/dp)s is the electrostrictive coupling 
constant, po is the unpertubed density of medium 
Consider the phase conjugation with a random phase 
plate which produces the chaotic transverse modulation 
of the phase of the incident E p with a characteristic size 
of the transverse inhomogeneity of the order of 5 — 50/j.m 
[l3l ]. In such a case the complex envelope amplitude of 
the inhomogeneous speckle field E p at the entrance to 
PC-mirror is given as a multimode random field [T(ij |: 

E p (r,z = 0)*E p ° A ^y 

0<j a ,j y <N g 

xexp[i2w{^K jx + ^-^k,. + i0 jx , jy }}, (3) 

Px Py 

where random phases Oj^j are the random numbers 
from interval [0, tt], are the real amplitudes of the 

spatial harmonics, p x ,Py are maximal transverse sizes (in 
X, Y-plane), f = (x,y) = (r,<j)), j x ,j y are integers cor- 
responding to N g — (16, 32, 64) plane waves with ran- 
dom phases 8j X j y . The random numbers Kj x = p x /*Pj X 
and Kj y = p v /pjy having the uniform distribution in the 
small vicinity of a 1 are responsible for a random tilt of a 
plane waves constituting the spatial Fourier spectrum of 
the light transmitted through a random phase plate. For 
the paraxial beam propagation the amplitude and phase 
structure of the complex field E p have the following form 
[r^ | in an arbitrary Z-plane : 

E p (f, z > 0) w E p ° exp (ikpz) 2J ex P [ * ^j x ,j v ] 

0<3x,jy<N g 

Aj „ tjy exp[i2K{^ + + ^(#2 + ^4)}].(4) 
Pjx Pjy ~ h 'f P jx P jy 

The interference pattern I(f, z) is responsible for the 
sound grating in a Brillouin PC-mirror. For the grating 
produced by a linear superposition of the random plane 
waves we have: 



FIG. 2: The plot of the light intensity maxima at a given 
moment i in a small volume inside the PC-mirror. The four X<m) 



pairs of spiral interference patterns with opposite handedness 
are located randomly in space. Their diameters and direc- 
tions are changed smoothly due to diffraction. The step of 
longitudinal modulation cos((k p + k a )z± 2<f>) is enlarged here 




I(f, z > 0) « \E p (r, z) + E s (r, z)\ 2 , (5) 
where the phase-conjugated Stockes wave E 8 is given by: 

E s (r, z > 0) « E p °exp (-ik s z) exp[—i jx , jy ] 

0<jx,jy<N g 

A jx , jy exp[-i2K {^ + y Jy+ -L(^4 + J4)}],(6) 

Pjx Pjy /' , , />,,, 

provided that E s w _E p * ( noiseless PC-mirror). Such an 
approach presumes the ideally perfect phase-conjugation 
with ultimate reflection R = 1 or interference of the non- 
interacting fields E s ~ E p * when the right hand sizes of 
Eq. (|1I2[) are equal to zero. 

The intensity snapshots given by ([5]) were obtained 
numerically by summation of the series (U) and © on 
standard Intel platform with dual — core 1.86 Ghz pro- 
cessor and 1 Gb memory using standard educational soft- 
ware. The numerical simulations show that optical vor- 
tex lines intertwine and form ropes. The average length 
of a rope is of the order of Rayleigh range D 2 /X. This 
may be interpreted as a consequence of diffractive di- 
vergence within the angle X/D (fig[T]). The longitudinal 
modulation by the Bragg factor cos[(k p + k s )z ± 2cj>] is 
accompanied by an additional bending and twisting at 
characteristic wavelength Xk- The latter is of the order 
of a several tenth of optical A, resembling Kelvin modes of 
the vortex lines [HI, UJ. The handedness of helical ropes 
changes randomly from one entity to another: the clock- 
wise and counterclockwise ropes appear with the equal 



FIG. 3: Gray scale intensity plot for the variational anzatz 
substituted in ([5]). The 128 x 2048 points arrays are 
used for the intensity distributions |S) in X x Z planes at 
fixed Y=64. The single sinusoidally modulated LGo 1 vortex 
line (see eq. ©) is shown for M = 0.1 ■ D (a) and M — D 
(b). Two LGo 1 distant vortex lines LGo 1 (c) with parallel (d) 
opposite (e) topological charges ±£ constitute the rope when 
the distance between them is reduced to zero(d,e). The step 
of longitudinal modulation is enlarged here 32 times. 

probability. Despite the statistical nature of a speckle 
pattern (figdj) the ropes are structurally stable: they 
appear at a different locations and have a variable sizes 
but their morphology is reproduced from a given statis- 
tical realization to another one. 

The physical interpretation of these numerically ob- 
served patterns is given in the figE] using the ([5]) with 
the electrical fields E p — E s * in the form of a two over- 
lapping elementary LGq 1 optical vortices: 

E p w exp (ikpZ - {xi 2 + yi 2 ) + Ufa) ■ y/ (xi 2 + yi 2 ) 
+exp (ikpZ - (x 2 2 + y2 2 ) ± Ufa) • \J (x 2 2 + y 2 2 ) , (7) 

where the location of the vortex cores (x^D — x + 
M ■ cos(2t:z/ Xk), {yiD = y + M ■ sin(2iTZ / Xk )) and 
(x 2 D = x + M ■ cos(2irz/X K + tt/2)), (y 2 D = y + M ■ 
sin{2-Kz I Xk + tt/2)) changes sinusoidally with period Xk 
along Z-eods, fa and fa are the local azimuthal angles 
around vortex cores. Such vortex-vortex (fig[3]i) and 
vortex-antivortex (fig(2t) pair interference pattern is used 
here as a variational anzatz in cq.([5]). The fig [3] repro- 
duces qualitatively a some features of the numerically 
generated speckle patterns (figHJ. 

In summary the nontrivial topology of the mul- 
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tiply connected optical speckle patterns was demon- 
strated numerically. Using the interference with the 
counter-propagating phase-conjugated speckle field we 
have shown the hidden twisted geometry of the multi- 
mode wave field composed of the randomly tilted plane 
waves fl6j l. Noteworthy the ropes exist without phase- 
conjugated counterpart: the reflection from PC-mirror 
makes these twisted entities visible due to the charac- 
teristic 7r-phase-jump in between the adjacent Bragg's 
sinusoidally -modulated roll patterns (figHJ). Thus it 
seems probable that a rope entities exist in a wide class 
of a superpositional physical fields like an electromag- 
netic or acoustical ones [Tl|. The other most evident 
examples are the blackbody radiation field in a cavity 
or a monochromatic field in a cavity with a rough mir- 
rors. Taking into account the phenomena of a so-called 
nonlinear superposition which take place for collisions 
of the optical solitons, vortices [2(| and formation of the 



optical vortex lattices [2l| one might expect the finding of 
a stable ropes in a nonlinear fields. The relevant exam- 
ples are photorefractive PC-mirrors 0, ultracold matter, 
e.g. phase-conjugating PC-mirrors based on nondegen- 
erate four- wave mixing in a trapped gases [!, Q and the 
situations of the slow [H, EJ [HJ and the fast light [23|, 
propagation in a resonant medium. 

The reported result extends the set of the possible 
forms of the photon's wavefunction. In addition to the 
well known photon wavefunctions in the form of the ele- 
mentary optical vortex lines, e.g. LG beams which pos- 
sess a helical wavefronts 0, H| , the twisted entities each 
composed of the several vortex lines offer a new form of 
the photon's wavefunction [25|, H|| . 
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